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ABSTRACT 

A parameterization of NAT (nitric acid trihydrate) 
clouds is developed for use in 2D models of the strato- 
sphere. The parameterization uses model distributions 
of HNC> 3 and H 2 O to determine critical temperatures for 
NAT formation as a function of latitude and pressure. Na- 
tional Meteorological Center temperature fields are then 
used to determine monthly temperature frequency distri- 
butions, also as a function of latitude and pressure. The 
fractions of these distributions which fall below the crit- 
ical temperatures for NAT formation are then used to 
determine the NAT cloud surface area density for each 
location in the model grid. By specifying heterogeneous 
reaction rates as functions of the surface area density, it 
is then possible to assess the effects of the NAT clouds 
on model constituent distributions. We also consider the 
increase in the NAT cloud formation in the presence of a 
fleet of stratospheric aircraft. The stratospheric aircraft 
NO x and H 2 O perturbations result in increased HNOjas 
well as H 2 O. This increases the probability of NAT for- 
mation substantially, especially if it is assumed that the 
aircraft perturbations are confined to a corridor region. 


1. INTRODUCTION 

Two dimensional photochemical models of the strato- 
sphere are important tools for the study of stratospheric 
dynamics and chemistry. They are also the major method 
for evaluating human influences 011 the stratosphere be- 
cause they require less computer time to run and contain 
a more complete photochemistry than the more compli- 
cated three dimensional models. In particular, they are 
being used to assess the potential effects of a fleet of 
stratospheric aircraft on stratospheric ozone [e.g., Jack- 
man et ai, 1991], It has been found that model predic- 
tions of aircraft-induced ozone depiction differ substan- 
tially, depending on whether or not the models include 
a parameterization of sulfate aerosol heterogeneous pro- 
cesses. Purely gas-phase models generally predict a de- 
pletion on the order of several percent,. If stratospheric 
sulfate aerosol effects are added, the predicted depletions 
are much less, or actually result in a slight increase in 
stratospheric ozone [Weiseristein et al. , 1991]. 


Heterogeneous reactions are also catalyzed by NAT 
clouds in the stratosphere. NAT clouds are different from 
the sulfate aerosols in that the temperature at which they 
will form depends on the amount of HNOsand H 2 0 in the 
stratosphere, and are not a ubiquitous feature of the lower 
stratosphere. NAT clouds are essential components of 
polar stratospheric chemistry and the properties of such 
clouds have been studied extensively [e.g., Turco et al . , 
1989]. 

Implementing NAT cloud distributions and the het- 
erogeneous reactions they catalyze into zonally averaged 
models of the stratosphere is difficult. This is becuase 
zonally averaged temperatures are generally too high to 
allow for NAT clouds to form. If one denotes by (c) the 
zonally averaged cloud surface area density, and the zon- 
ally averaged temperature by (X), then (c(X)) ^ c({T}). 
To get around this difficulty while remaining consistent 
with the philosopy of a 2D model, we use temperature 
probability distributions characterizing a particular lati- 
tude and pressure level in the model. These distributions 
are obtained using National Meteorological Center tem- 
perature data. 

NAT cloud formation depends on the concentrations 
of HNOijand H 2 O in addition to the temperature. Be- 
cause stratospheric aircraft inject both into the strato- 
sphere, it is likely that a stratospheric aircraft fleet will 
increase the probability of NAT cloud formation. Peter et 
al. [1991] attempted to assess the magnitude of the effect 
at 70°N, and concluded that the probability of finding an 
NAT cloud at that latitude might be doubled with the 
introduction of such a fleet. This paper is in the same 
spirit as their work, but attempts to develop an annual 
NAT surface area density climatology for the latitudes 
and pressure levels of the GSFC 2D model. We do this 
both with and without the assumption of a stratospheric 
aircraft fleet in order to allow a 2D model assessment of 
the chemical effects of stratospheric aircraft on the atmo- 
sphere. 

2. CRITICAL TEMPERATURES FOR NAT FORMA- 
TION 

According to Hanson and Mauersherger. [19SS], the 
relationship between the NAT saturation temperature. 
X S1 the partial pressure of HNO ;1 . Pus o ;( . and the par- 
tial pressure of H^O. P\\ 7 o. is given by. 
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log 10 ( p HNO3) = m(T a ) log lfl (P H 2o) + KTs) 


( 1 ) 


and 

m{T,) = -2.7836 - 0.00088T, (2) 

1 1 3Q7 0 

b(T,) = 38.9855 + °- 009179T ^ ( 3 ) 

The stated range of validity of this equation is for tem- 
peratures between 180 K and 200 K. In this study, we 
assume that the equation is valid up to 210 K, and that 
NAT clouds will form wherever the temperature is below 
T s . This assumption may not be correct, however. Peter 
et a/.[1991] have discussed the possibility that a substan- 
tial supersaturation on the order of 3 K is required to 
nucleate the NAT clouds. If this is the case, then using 
Eq. 1 will result in a prediction of the critical tempera- 
tures for NAT cloud condensation that is too high. 

Given Eq. 1, and model distributions of HN0 3 and 
H 2 0, a critical temperature distribution can be calculated 
as a function of latitude and pressure, as shown in Fig- 
ure 1 . In this case, a subsonic aircraft fleet has been as- 
sumed, but no stratospheric aircraft fleet. The addition of 
stratospheric aircraft increases the critical temperatures. 
The amount of increase depends on the assumptions made 
about the zonal distribution of HNO 3 . For instance, the 
dynamics of the stratosphere may be such that aircraft 
corridors are formed inside of which the consituent con- 
centrations are significantly greater than the zonal mean 
values. Tracer simulations using the GSFC 3D chemistry 
and transport model and European Center for Medium 
Range Weather Forcasting (ECMWF) wind fields suggest 
maximum concentrations of up to 5 times the zonal mean 
during the winter [Douglass et < 1 /., 1992). If one considers 
all of the injected HN0 3 and H 2 0 from the stratospheric 
aircraft to be concentrated in a corridor with a width 
of 72° longitude, the local concentration increase inside 
the corridor will be 5 times the perturbation to the zonal 
mean, while outside the corridor the contribution from 
the aircraft fleet will be 0 . Using this assumption, the 
increase in NAT formation temperatures can be up to 6 
I\, as shown in Figure 2. 

3. TEMPERATURE PROBABILITY DISTRIBUTIONS 

As mentioned above, zonally- averaged temperatures 
are generally too warm to dip below the NAT critical tem- 
peratures. if one considers the longitudinal distribution 
of temperatures at a particular latitude, however, it could 
be that a substantial fraction of the band is below the con- 
densation temperature, and could support NAT cloud for- 
mation. To quantify this effect, we took a month s worth 
of NMC temperature data, and binned the temperatures 
into lOMatitudo bands. At each NMC pressure level, the 
temperatures in the band were used to produce a tem- 
perature probability density distribution for each month, 
Pn(T 0, ./>,), where 6 t = ( -85, -75, 85) are the lat- 
itude coordinates and { Pj } are the NMC pressure lev- 
els. This distribution is normalized such that 1 (U gnes 
the probability of finding a temperature between T and 
7 + (IT. For tiiis paper, NMC temperatures for one year, 
from Mav. 1991 to April, 1992, were used. By consider- 



Fig. 1 Critical Temperatures for NAT formation: As- 
suming the relationship between saturation tem- 
perature and distributions of HN0 3 and H 2 0 
given in Hanson and Mauersberger [1988], criti- 
cal temperature distributions can be calculated. 
Shown is the January, 1992 distribution, where 
the HN0 3 and H 2 0 fields were taken from a 
model run including sulfate aerosol heterogenous 
processes and a subsonic aircraft fleet. 



Fig. 2 Increase in critical temperatures from addition 
of a stratospheric aircraft fleet: The increases in 
stratospheric HN0 3 and H 2 0 from stratospheric 
aircraft can substantially increase the critical 
temperatures for NAT formation. This effect is 
on the order of 2 K if the aircraft perturbation 
is assumed to be distributed zonally. If the air- 
craft perturbation is assumed to be confined to a 
tight corridor, concentrations are higher and so 
is the temperature. The figure above was calcu- 
lated assuming an aircraft corridor width of 72° 
longitude. 
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ing more than one year of NMC data, a better estimate 
of the climatological temperature distributions could be 
obtained. 

With the critical temperatures calculated using the 
GSFC 2D model and the NMC temperature distributions, 
the probability that NAT clouds will be found can be 
calculated as a function of month, latitude, and pressure. 
Denoting this fraction by /, 

fm{T c ,8i,pj) = / Pm(T, 8i,pj)dT, (4) 

Jo 

where T c is the condensation temperature for NAT clouds 
(which we assume is equal to the Hanson and Mauers- 
berger saturation temperature, Eqs. 1-3, above). As 
shown in Figure 3, which depicts the calculation of / at 50 
mb, 65° for January, 1992, / can be a rapidly increasing 
function of T c . If T c increases slightly due to the introduc- 
tion of stratospheric aircraft, / can increase substantially. 

A plot of / for January 1992 is shown in Figure 4. 
In the case shown, no stratospheric aircraft fleet has been 
assumed. For this month, the peak probability that an 
NAT cloud will form occurs in the polar region, and has 
a magnitude of .2. From the figure, we might expect to 
see an NAT cloud forming at, say, 70° and 50 mb with 
a probability of about 0.1, or roughly three days out of 
the month. This appears to be in fair agreement with 
the results of the study of Peter et al. [1991], which also 
calculated such probabilities. The figure also indicates 
a substantial probability that NAT clouds will form in 



Fig. 3 Calculation of / from T r and P(T): The solid 
curve at the bottom of the figure is the tem- 
perature probability distribution at 65° N and 50 
mb. the dashed line is the integral of this dis- 
tribution up to the temperature T. The vertical 
dashed lines show the critical temperatures for 
NAT condensation inside and outside of an air- 
craft corridor with a longitudinal width of 72 (> . 
Because / can increase rapidly, f(T r ) inside the 
corridor can be substantially larger than f(T r ) 
outside of the corridor. 
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Fig. 4 January 1992 NAT formation probabilities: No 
stratospheric aircraft are included in this cal- 
culation. The peak in the January formation 
probability occurs near the north pole, at ap- 
proximately 70 mb. A significant tropical for- 
mation probability is also seen, which peaks a 
bit lower down at approximately 100 mb. Very 
small probabilities of NAT formation are found 
in the northern midlatitudes. Temperatures in 
the southern hemisphere are too warm for NAT 
clouds to exist there. 
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Fig. 5 Increase in January NAT formation probability 
from the addition of a stratospheric aircraft fleet. 
Shown is the increase assuming a 72° aircraft 
corridor, inside the corridor. 


the tropics, peaking at about 100 mb, a little lower than 
the polar peak. The probability of cloud occurence is a 
minimum at northern midlatitudes, and is absent in the 
summer hemisphere, reflecting the warmer temperatures 
there. 

When the N() x and HjO perturbations from a strato- 
spheric aircraft fleet are included in the calculations, the 
resulting probabilities of NAT cloud occurrence are in- 
creased. The increase in probability of cloud formation 
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Fig. 6 Seasonal variation of NAT formation probabil- 
ity: Shown is the' variation of / as a function 
of day and latitude, on the 90 mb surface. No 
stratospheric aircraft perturbation is included. 
The largest probabilities are in the south po- 
lar regions in the months of July, August, and 
September. Substantial probabilities of occur- 
rence are also found in December/ January in 
the northern high latitudes, and in the tropics 
in February/ March. At 50 mb, iio tropical cloud 
formation is predicted, but the north and south 
polar distributions are similar to those depicted 
here. 

has somewhat the same morphology as / itself, and has 
a similar magnitude. If the perturbations are assumed to 
be confined to a corridor, the increases are larger than if 
the perturbations are assumed to be distributed zonally 
(Figure 5). 

Distributions similar to Figure 4 can be obtained 
for each month. Figure 6 depicts the seasonal variation 
of NAT cloud occurrence probability on the 90 mb sur- 
face. The north polar peak in occurrence probability oc- 
curs in December/ January. The south polar occurrence 
probability peaks in August, and is large in both July 
and September. Finally, the tropical probability peaks 
in February/ March, and is at a minimum in October. If 
one moves to the 50 mb surface, one finds only very small 
probabilities of tropical cloud formation, while the north- 
ern and southern polar regions remain about the same. 

If we assume that the NAT cloud that forms has a 
value of 10“ 8 cm -1 , as given by Turco et al, [1989], the 
distributions define a seasonally- varying NAT cloud sur- 
face area density distribution. (This assumption may in- 
troduce a substantial uncertainty into the calculation. It 
should be adequate for a 2D sensitivity study.) The sur- 
face area density distribution can then be introduced into 
the model, and used to calculate the effects of hetero- 
geneous reactions catalyzed by the NAT clouds on con- 
stituent distributions. Since the NAT cloud freqency is 
increased with the addition of a stratospheric aircraft 
fleet, it appears that including NAT cloud processing in 
assessments of potential stratospheric aircraft effects on 


the stratosphere may result in important changes. These 
changes remain to be assessed. 
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